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Abstract 

Einstein's vacuum equations are solved up to the second approximation 
for imploding quadrupole gravitational waves. The implosion generates a 
black hole singularity irrespective of the strength of the waves. 



1 Introduction 

Abraham and Evans [1] [2] have reported numerical calculations describing the 
implosion of quadrupole gravitational waves. They find that the waves disperse 
if sufficiently weak, but form a singularity if strong enough. This work has im- 
plications for the formation of black holes and for critical phenomena in general 
relativity. 

In this paper we describe analytical approximate solutions for incoming 
quadrupole waves which lead to a different conclusion. We first give the solution 
of the linear approximation to Einstein's equations corresponding to an axially 
symmetric incoming quadrupole wave. Then, using previously published work, 
we can write down the second approximation exactly. We find that there is a 
singularity after the implosion of the wave, whatever its strength. 

Most work on gravitational waves in the past has, naturally, dealt with 
outgoing waves, and quadrupole waves have often been studied. For reviews see 
[3] [4] [5]. In this paper we shall use the double-series approximation method 
[6] [7] which was set up for outgoing waves, but can be adapted to deal with 
incoming waves. We are therefore able to use previous work, particularly [7], 
with hardly more than a few changes of sign, to get the results of this paper. 

In Section 2 we introduce the metric we use to study the imploding 
waves and in Section 3 we describe our approximation method. The linear and 
second approximations are given in Sections 4 and 5, and there is a concluding 
Section 6. 



1 



2 The metric 

We confine ourselves to the axially symmetric case, and use a metric of the form 
given by Bondi [8] [9] . For the study of outgoing waves one can take 

ds 2 = -r 2 (bd9 2 + csin 2 8d<j) 2 ) + edu 2 + 2fdrdu + 2rgd9du, (1) 

where c = b^ 1 and b, e, /, g are functions of r, 9 and u. The flat form of this is 

ds 2 = -r 2 (d9 2 + sin 2 6d(f) 2 ) + du 2 + 2drdu (2) 

from which we see that u is a retarded time, obtained from the t of the Minkowski 
metric by putting u = t — r. The Schwarzschild solution in these coordinates is 

ds 2 = -r 2 (d9 2 + sin 2 6d<j) 2 ) + (1 - Imr-^du 2 + 2drdu. (3) 

In this paper we are concerned with incoming waves. We therefore use 
an advanced time v given in terms of the ordinary Minkowski t by v = t + r. 
The flat and Schwarzschild metrics now become 

ds 2 = -r 2 (d9 2 + sin 2 9d(f) 2 ) + dv 2 - 2drdv, (4) 
ds 2 = -r 2 (d9 2 + sin 2 9d<j) 2 ) + (1 - 2mr- 1 )dv 2 - 2drdv. (5) 

which can be obtained from (2) and (3) by the coordinate transformation 

u = —v, (6) 

r, 9, 4> being unchanged. (For a discussion of the Schwarzschild solution in ad- 
vanced and retarded time coordinates see [10]). 

In what follows we shall use (6) to transform known solutions for out- 
going waves to new solutions for incoming waves. The metric for the incoming 
waves will be written 

ds 2 = -r 2 (Bd9 2 + Csin 2 9d<p 2 ) + Ddv 2 + 2Fdrdv + 2rGd9dv, (7) 

where C = B' 1 and B,D,F,G are functions of r,9,v. Metrics (4) and (5) are 
special cases of (7). 

3 The approximation method 

In the double-series approximation method as used heretofore one seeks solu- 
tions of Einstein's vacuum equations for an isolated source by expanding the gik 
in terms of two parameters m and a, referring respectively to the source's mass 
and a characteristic length associated with the source: 

9ik =J2J2 mPaS *k, (8) 
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where 9 ik are independent of m and a. The flat background metric is obtained 
by putting p = 0, s = 0, the Schwarzschild metric in Bondi coordinates comes 
from p = 1, s = 0, and the linear approximation consists of the infinite series 
of terms p = 1, s > 0. The latter comprise the multipole terms, and these can 
be prescribed arbitrarily. By considering non-linear terms p > 2, s > 1 one 
can study the interactions of the multipole effects with themselves, or with the 
Schwarzschild mass term. For further details, see [6]. The linear quadrupole 
terms are those with p = 1, s = 2. 

For incoming waves as considered in this paper the concept of an iso- 
lated source is not applicable because the waves are coming from infinity. How- 
ever, we can still consider the radiation in terms of multipolcs, and study their 
interactions. In this paper we are interested only in quadrupole waves, and 
their interactions with themselves. There is no (10) (Schwarzschild) term, and 
no other multipole terms. We can therefore simplify the notation by putting 
A = ma 2 , and expanding the field equations in the single parameter A. Thus 
(8) can be written 

9 ik = 22 Xn ( 9 ) 

n=0 

(n) 

where the 9 ik are independent of A. For the metric (7) we have 

°° in) 

- 922 =r 2 B = r 2 [l + ]TA" B ], (10) 

n=l 

°° («) 

-533 = r 2 Csm 2 9 = r 2 sin 2 9[l + ^ A" C], (11) 

71=1 
°° (™) 

944 = D = l + ]TA n D, (12) 

n=l 

oo , , 

9i4 = F = -1 + 2 A" F , (13) 

71=1 

OO , . 

924 = rG = r^rG, (14) 

71=1 

In this paper we confine ourselves to the terms with n = 1 and n = 2. 



4 The linear approximation 

We first substitute (7) into the vacuum equations 

R ik = (15) 
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and separate the linear approximation which, because of the axial symmetry, 
consists of seven second order equations 

*lm(flifc)=0. (16) 

The solution of these for incoming quadrupole waves can be written down from 
the corresponding solution for outgoing waves given in [6] equations (6.6), (6.9), 
(6.10) and (6.11), by the use of transformation (6). It is 

D = -P 2 (2r- 1 h-2r- 2 h + r- 3 h) 1 (17) 

W , ^ 

F = 0, (18) 
G = (l/6)P2(-2r- 1 h-4r- 2 h + 3r- 3 h), (19) 
jB = (1/2) sin 2 6(r- v h + r- 3 h), (20) 

<8 - -g, (21) 

where P 2 denotes the second Legendre polynomial, / denotes differentiation 
with respect to 9, Xh(v) is the quadrupole moment, and an overhead dot means 
differentiation with respect to v. We shall assume that 

h = 0,v<vi; h = 0,v>v 2 (22) 

and that h is a smooth function of v in vi < v < v 2 . Thus the incoming 
quadrupole wave is of finite duration. The spacetime is flat for v < v\ . In the 
linear approximation (17)- (21) it is flat also for v > v 2 , but this is not so for the 
second approximation, as we shall see. 

The functions D, G, B are singular at r — 0, and one can confirm that 
r = is a curvature singularity. It therefore appears that in the linear ap- 
proximation the incoming quadrupole waves generate a physical singularity in 
the neighbourhood of r = 0, whatever their strength. However, this singular- 
ity disappears at v — v 2 because of (22). This non-physical behaviour will be 
corrected in the second approximation. 

5 The second approximation 

We now have to solve seven equations of the form 

(2) (1) 

<S>im{9 tk ) = ^i m (9 lk ), (23) 

(2) 

where the left-hand side is linear in 9 ik (and their derivatives). The right-hand 

(i) 

side is non-linear in the 9 ik (and their derivatives), which are known from the 
linear approximation. 
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This problem was completely solved for outgoing quadrupole waves in 
[7] and all we need to do to obtain a solution for incoming quadrupole waves is 
to apply the transformation (6) to that work. The result of doing this is given in 
the Appendix. Here we omit transient terms, i.e. those which vanish for v > 
because our interest is in the formation of permanent singularities. The result 
is 

(2) 1 

B = -^ sin M 2 + sin 2 fl)y, (24) 

(2) (2) 

C = - B, (25) 

(2) 1 

D = w r, (26) 

< 2 ) , ^ 

F = 0, (27) 

(2) 1 

G = — sin (9 cos (9(4 + 3 sin 2 <9)Y", (28) 

where s = sin#, c = cos 9, and 

h dv. (29) 

Y(v) in (24)-(28) vanishes for v < v\ but is a negative constant for 
v > V2- This represents a permanent change in the metric. In fact, by a 
coordinate transformation similar to (9.13) of [6], one can show that the metric 
for v > V2 is, up to order A 2 , that of Schwarzschild for a particle of mass 
A 2 | Y(v 2 ) | /30. The result is valid whatever the value of A, which means that 
it applies whether the waves are weak or strong. 



6 Conclusion 

We have shown that an incoming quadrupole wave produces, in the second 
approximation, a singularity which persists after the wave has ceased. This 
result is independent of the strength of the wave. At this level of approximation 
the singularity is a black hole, since the metric has Schwarzschild form. 

The solution we have presented here is, so to speak, the mirror image 
of the corresponding terms of the solution for outgoing waves. In [6] and [7] Y 

(24) 

(there denoted by Y ) traces a loss of mass of an isolated source, represented 
by a Schwarzschild particle ; this corresponds to the energy carried away by the 
waves. In the imploding case a corresponding amount of energy flows inward, 
converging on r = 0. 

Our result, that a black hole is formed whatever the strength of the 
incoming wave, disagrees with that of [1] and [2] . One obvious possible cause of 
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the disagreement is that whereas our work stops at the second approximation 
the numerical calculations of Abraham and Evans include contributions from 
higher approximations. However, this explanation seems implausible for the 
following reason. The disagreement occurs when the incoming wave is weak, i.e. 
for small values of A; but for these values one would expect convergence to be 
quicker, and the higher powers in the expansion to be less important. 



We gratefully acknowledge an interesting correspondence with Dr. John 
M Stewart. 



Appendix 

Here we give the complete solution to the second approximation summarised in 
Section 5. It is obtained by using the transformation (6) on the corresponding 
solution of the second approximation for the outgoing wave problem, namely 
(3.7)-(3.11) of [7]. The only matter which is not quite straightforward is the 

24 

relation of Y of [7] and Y in Section 5. The essential point here is that the 
former is defined by 

(24) 

= {ti") 2 (30) 



d Y , «... ■> 



du 

so on applying the transformation we have 



We are now at liberty to integrate this equation: 

r ...2 

y = - h dv (32) 

J — oo 

the lower limit being chosen because spacetimc is assumed to be empty and 

(24) _ 

static at v — — oo. Thus Y is replaced by Y defined in (29). 
The complete second approximation is 

(2) 

D = [(§ - 3s 2 + f s 4 )hh" + (f - 9s 2 + f s 4 )hh +±Y]r- 1 
+ [(-2 + 10s 2 - f s 4 )hh + (-2 + 9.s 2 - f s 4 )h 2 }r- 2 

+ (! - f * 2 + f^ 4 )^ 3 + [(I- 2 - %°*)hh 

+ X s ' - ^§s 4 )h 2 ]r- 4 + (§- 6s 2 + %s 4 )hhr- 5 



6 



+(i _ 3s 2 + f S 4 )^ 2 r- 6 (33) 

(2) 

p = s 4 [^h 2 r- 2 + ^hhr~ 4 + ^h 2 r- 6 } (34) 
G = [(-±cs + %cs 3 )hh™ + (-lcs+§cs 3 )hh -(^cs + ^cs^r- 1 
+[(-2cs + lcs 3 )hti + {-2cs + f cs 3 )h 2 ]r- 2 

+ (^ cs _ M.cs 3 )hhr- 5 + (fcs - f cs 3 )h 2 r- e (35) 
§ = [(-l s 2 + ^)hh" + (-ls 2 + l oS 4 )hh +(_^__i_ a 4 )y]r -i 

+ 1^-2 + ( _9 S 2 + ^ s 4)^ r -3 + [1 ^ + ( 75 g 2 _ 175^2^-4 
+ (-f S 2 + 69^)^-5 + (_7 S 2 + W s 4j fe 2 r -6 (gg) 

(2) 

C = K^ 2 -4lJ s4 )^ IV + (f s2 -|j s4 )^ +(M s2 +i3o s4 )^ _1 

+ 1^-2 + ( 9 S 2 _ M a 4)^ r -3 + [1,4^ + ( _ 7| s 2 + 175 ^2^-4 
+ (^ s 2 _ 69 s 4)^ r -5 + ( 7 S 2 _ 1^4^-6 (3?) 

where s = sin and c = cos 0. 
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